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Abstract 

We study the reduced Bogoliubov-Dirac-Fock (BDF) energy which al- 
lows to describe relativistic electrons interacting with the Dirac sea, in 
an external electrostatic potential. The model can be seen as a mean- 
field approximation of Quantum Electrodynamics (QED) where photons 
and the so-called exchange term are neglected. A state of the system is 
described by its one-body density matrix, an infinite rank self-adjoint op- 
erator which is a compact perturbation of the negative spectral projector 
of the free Dirac operator (the Dirac sea). 

We study the minimization of the reduced BDF energy under a charge 
constraint. We prove the existence of minimizers for a large range of values 
of the charge, and any positive value of the coupling constant a. Our 
result covers neutral and positively charged molecules, provided that the 
positive charge is not large enough to create electron-positron pairs. We 
also prove that the density of any minimizer is an L 1 function and compute 
the effective charge of the system, recovering the usual renormalization of 
charge: the physical coupling constant is related to a by the formula 
«ph ys — a(l + 2«/(37r) log A) -1 , where A is the ultraviolet cut-off. We 
eventually prove an estimate on the highest number of electrons which 
can be bound by a nucleus of charge Z. In the nonrelativistic limit, we 
obtain that this number is < 2Z, recovering a result of Lieb. 

This work is based on a series of papers by Hainzl, Lewin, Sere and 
Solovej on the mean-field approximation of no-photon QED. 



1 Introduction 

In this paper, we study a model of Quantum Electrodynamics (QED) allowing to 
describe the behavior of relativistic electrons in an external field and interacting 
with the virtual electrons of the Dirac sea, in a mean-field type theory. This work 
should be seen as the continuation of previous papers by Hainzl, Lewin, Sere and 
Solovej [T^] — [IB], in which a more complicated model called Bogoliubov-Dirac- 
Fock (BDF) is considered. This project was mainly inspired of an important 
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physical paper by Chaix and Iracane [51 [5] in which a model of the same kind 
was first proposed. We start by summarizing the physical motivation before 
defining the model properly. 

Dirac introduced his operator in 1928 [7] with the purpose to describe the 
behavior of relativistic electrons. It is defined as 

3 

D° = a kd k + 13 := -ia ■ V + (3 (1) 

k=l 

where a — (a\, ct2, a^) and (3 are the 4x4 Dirac matrices [27]. The operator D° 
acts on L 2 (R 3 , C 4 ). Contrary to the non-relativistic Hamiltonian —A/2, the op- 
erator D is unbounded from below: a(D°) = (— oo, — 1] U [1, oo). This property 
is known to be the basic explanation of various peculiar physical phenomena 
like the possible creation of electron-positron pairs or the polarization of the 
vacuum. The model that we shall study is a rough approximation of Quantum 
Electrodynamics but it is able to reproduce many of these physical phenomena. 
We refer to [T2] — [TS] for more details. 

In QED, one can write a formal Hamiltonian acting on the usual fermionic 
Fock space, in Coulomb gauge and neglecting photons [E3 Eq. (1)]. The mean- 
field approximation then consists in restricting formally this Hamiltonian to a 
special subclass of states in the Fock space, called the Hartree-Fock states. Any 
of these states is uniquely determined by its one-body density matrix which is a 
self-adjoint operator < P < 1 acting on L 2 (R 3 , C 4 ). Often P is an orthogonal 
projector. The QED energy then becomes a nonlinear functional in the variable 
P, which can be formally written as follows 



-QED 



(P) = tr(P°(P - 1/2)) - a [( V{x)P ^ {v) d X d y 

P[P -^) P[P - m {y) dxdy ^ ff l(P-i/2)(^)|^ 



2 JJwxW \ x ~ V\ ' 2 JJmzxRz \x-y\ 

(2) 

where for any operator Q acting on L 2 (M 3 ,C 4 ) with kernel Q(x,y), pq is for- 
mally defined as pq(x) = trci(Q(x,x)). Recall Q(x,y) acts on 4-spinors, i.e. is 
a 4 x 4 complex hermitian matrix. The first term of @ is the kinetic energy of 
the particles, whereas the second term describes the interaction with an external 
electrostatic field created by a smooth distribution of charge v (describing for 
instance a system of classical nuclei). The last two terms account for the inter- 
action between the particles themselves. We have chosen a system of units such 
that h = c = 1, and also such that the mass m e of the electron is normalized to 
1. The constant a — e 2 (where e is the bare charge of an electron) is a small 
number called the Sommerfeld fine- structure constant. 

Expression ^ is purely formal: when P is an orthogonal projector on 
L 2 (M 3 , C 4 ), P — 1/2 is never compact and none of the terms above makes sense 
a priori. However, it is possible to give a meaning to ([2|) by restricting the 
system to a box and imposing an ultraviolet cut-off. One can then study the 
thermodynamic limit, i.e. the behavior of the energy and of the minimizers 
when the size of the box goes to infinity (but the ultraviolet cut-off is fixed). 
This approach was the main purpose of [15j . 
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The last two terms of © are respectively called the direct term and the 
exchange term. In theoretical studies of the Hartree-Fock model, the exchange 
term is sometimes neglected |26| . The above energy then becomes (formally) 
convex, a very interesting simplification both from a theoretical and numerical 
point of view. Refined models exist: in relativistic density functional theory 
for instance, the exchange term is approximated by a function of the density 
9{P-\li\ an d its derivatives only, see, e.g., the review [TT]. Neglecting the last 
term, one is led to consider the following reduced formal functional 

£ r - QED (P) = HD\P - 1/2)) -all V{x)P ^f y) d X dy 

+ a[f PiP-i/^pjp i,^y) 
ZJJrixr 3 \x-y\ 
As usual, one is interested in finding states having lowest energy, possibly in 
a specific subclass. In QED, a global minimizcr in the Fock space is interpreted 
as being the vacuum, whereas other states (containing a finite number q of real 
electrons for example) are obtained by assuming a charge constraint. When 
the external field vanishes (y = 0) and for any values of the coupling constant 
a > 0, one easily proves that £ °_q ED has a unique global minimizer which is the 
negative spectral projector of the free Dirac operator: 

P- ■= X(-oo.o](D°). 

The precise mathematical statement is that when the system is restricted to a 
box of size L with an ultraviolet cut-off A, the above energy is well-defined; it 
has a unique minimizer 

Pl = X(-oo,0]( D l) 

where P)\ is the Dirac operator acting on the box with periodic boundary con- 
ditions. The sequence P° converges (in a weak sense) to P° which is thus 
interpreted as the unique global minimizer of P i— > £®-Qet>(P)- ^ tne exchange 
term is not neglected, the situation is more complicated and we refer to |15j 
where the thermodynamic limit was carried out. 

The fact that P° is found to be the global minimizer of our formal energy is 
not physically surprising. This corresponds to the usual Dirac picture [7|,[8t[9t[T0] 
which consists in assuming that the vacuum should be seen as an infinite system 
of virtual particles occupying all the negative energy states of the free Dirac 
operator. Notice however that when the exchange term is taken into account, 
this picture is no longer valid: P° does not describe the free vacuum which is 
instead solution of a complicated translation-invariant nonlinear equation, see 

We want to emphasize the importance of the subtraction of half the identity 
in all the terms of the above energy ^ . Indeed, the kernel of the translation- 
invariant operator P_ — 1/2 is 

-3/2 f( „ „.\ „v. $( U \ _ D (k) 



(PH. - l/2)(x, y) = (2ir)-^ 2 f(x - y) where f(k) = - 



2\D°(k)l 



If we assume that there is a cut-off A in the Fourier domain, i.e. supp(/) C 
P(0, A), it is then possible to compute the density 

p [P o_ l/2] = (2^)- 3 / 2 tr C 4(/(0)) = (27T)- 3 I tr C 4(/>))dfc = 0, (4) 

JB(0,A) 
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the Dirac matrices being trace-less. We therefore obtain that the free vacuum 
has no density of charge, which is comforting physically. 

When the external field does not vanish, the main idea is then to subtract 
the (infinite) energy of the free vacuum S °_qed (P- ) to f5|), in order to obtain a 
finite quantity. This yields the so-called (formal) reduced-Bogoliubov-Dirac-Fock 
energy (rBDF) which was already studied in [13] and is more easily expressed 
in terms of the difference Q = P — P°_, 

£»(P-P°_) = «£» Q ed(P)-C Q ed(P-T 

= tr { D\P-P«))-a[[ U{X)p[P - P °- ]{y) dxd y 



\x - y\ 

p [P _ p a ] (x)p [P _pa ] (y) 



2 I a; - y\ 



dxdy. (5) 



Note that we have used (J4j) . What we have gained is that Q = P — P°_ can now 
be a compact operator (it will indeed be Hilbert-Schmidt). We recall that P is 
the density matrix of our Hartree-Fock state, hence it satisfies < P < 1 which 
translates on Q as -P° < Q < 1 - P° := P£. 

A (formal) global minimizer Q of is interpreted as the polarized vacuum 
in the presence of the external density v. Formally, it solves the self-consistent 
equation 

Q = X(-oo,o){Dq) - P- (f ~ 
D Q = D° + a{ PQ -v)*\-\-\ [b) 

In order to describe a physical system containing a finite number q of real 
electrons, it is necessary to minimize the above energy not on the full class of 
states, but rather in a chosen charge sector, i.e. over states satisfying the formal 
charge constraint "tr(Q) = tr(P — P°_) = q" . Then a minimizer will satisfy the 
following equation 



'<> 



Do = D Q + a(po-v)*\-\- 1 (7) 



where p is a Lagrange multiplier due to the charge constraint and interpreted 
as a chemical potential. The operator 5 is a finite rank operator satisfying 
< 5 < 1 and Ran(£) C ker(Z?Q — fi). Notice the number q does not need to 
be an integer as one may want to describe mixed states (in which case 6 ^ 0). 

We see that in both cases (minimization with or without a charge constraint), 
a minimizer always corresponds to filling energies of an effective Dirac operator 
up to some Fermi level p. This corresponds to original ideas of Dirac. For the 
general BDF theory, the idea that one can have a bounded below functional 
whose minimizer satisfies this kind of equation was first proposed by Chaix and 
Iracane OH]. 

In this paper, we shall prove that the range of g's such that minimizers exist 
is an interval [q m ,Q_M] C M which contains both the charge of the polarized 
vacuum (the global minimizer of the energy, solution of ©) denoted by qo, and 
Z = L3 v. This proves the existence of neutral molecules and of positively 
charged molecules the charge of which is not too big, because in this case one 
has qo — 0. This extends previous results proved for the BDF theory with 
the exchange term in [14] : sufficient conditions were given for the existence of 
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minimizcrs, but these conditions could only be checked in the nonrelativistic or 
the weak coupling limits. In the present paper, we shall also give interesting 
properties of a minimizer when it exists, and provide a bound on the maximal 
number of electrons which can be bound by a nucleus of charge Z, following 
ideas of Lieb [20] . 

The mathematical formulation and the proofs of the above statements are 
not straightforward. 

The first (and main) difficulty is that we do not expect that a solution Q 
of Equations ^ or (J7|) is a trace-class operator. Indeed our results below will 
imply that in most cases it cannot be trace- class. This is a big problem as in the 
energy ^ the first term is expressed as a trace, as well as the total charge of 
the system which we formally wrote "tr(Q)" in the previous paragraphs. This 
issue was solved in |12) where it was proposed to generalize the trace functional 
and to define the trace counted relatively to the free vacuum P° as 



As we shall see, any minimizer Q will have a finite so-defined P° -trace, which 
does not mean that Q is trace-class. 

If we do not expect Q to be trace-class, there is a problem in defining the den- 
sity of charge pq. Indeed it is known that in QED there are several divergences 
which need to be removed by means of an ultraviolet cut-off. In previous works 
[T2"] — [T§], a sharp cut-off A was imposed: the space L 2 (K 3 , C 4 ) was replaced by 
its subspace consisting of functions that have a Fourier transform with support 
in the ball of radius A. This allowed to give a solid mathematical meaning to 
the energy ([5]) . In [T^l [T3] , it was proved that the energy has a global minimizer 
Q, solution of In [13] , sufficient conditions were given on q to ensure the 
existence of a ground state in the charge sector q with the exchange term. They 
could only be checked in the nonrelativistic or the weak coupling limit. 

In this paper, we propose other kinds of cut-offs which seem better for ob- 
taining decay properties of the density of chargeQ. Essentially, they consist in 
replacing the Dirac operator D° by D^(p) = (a ■ p + /3)(1 + ((\p\ 2 / A 2 )) where 
£ is a smooth function growing fast enough at infinity. We call these cut-offs 
smooth in contrast to the previous sharp cut-off. But many of our results will 
also be valid in the sharp cut-off case. 

Even with an ultraviolet cut-off, a minimizer Q will in general not be trace- 
class. But we shall be able to prove that anyway its density of charge is an L 1 
function: pq 6 L 1 (R 3 ) . This information can then be used to prove the existence 
of all atoms and molecules which are either neutral or positively charged and 
do not have a too strong positive nuclear density. Also we shall prove a formula 
which relates the integral of pq and q — tr P o (Q) of the form 



(see Theorem U] for a precise statement depending on the chosen cut-off A) . 
When q Z, this proves that / ffi3 pq ^ q = tr P o(Q), hence Q cannot be 
trace-class. 

1 A similar remark was made in 12 1 1 in the context of non-relativistic QED. 



tr>o (Q) := tr(P°QP°) + tr(P° 



QP -). 




(8) 



5 



The fact that a minimizer is not trace-class but its density is anyway an L 1 
function can first be thought of as being a technical issue. But Equation © 
has a relevant physical interpretation. It means that the total observed charge 
/r3 PQ — Z is different from the real charge q — Z of the system. Hence the 
mathematical property that a minimizer is not trace-class is well interpreted 
physically in terms of charge renormalization. We even recover a standard 
charge renormalization formula in QED, see [HO Eq. (8)] and [TTj, Eq. (7.18)], 
although we use a simple model without photons and within the Hartree-Fock 
approximation with the exchange term removed. 

As announced before, we shall prove in this paper that minimizers exist if 
and only if q £ [SmiSAfL an interval which contains both Z and the charge go 
of the polarized vacuum. We shall also derive some bounds on q m and qm, 
assuming that the nuclear charge distribution is not too strong. Essentially we 
prove that q m < is very small and that 

Z < q M <2Z+ o (1). 

In the nonrelativistic limit we recover the usual bound of the reduced Hartree- 
Fock model which can be obtained by a method of Lieb [20] • 

In the next section, we define the reduced BDF energy (|5|) properly and state 
our main results. Proofs are given in Section [3J 

Acknowledgment. M.L. and E.S. acknowledge support from the ANR project 
"ACCQUAREL" of the French ministry of research. 

2 Model and main results 

In the whole paper, we denote by & p (Sj) the usual Schatten class of operators 
Q acting on a Hilbert space Sj and such that tr(|Q| p ) < oo. We use the notation 
Q ee := P®QP®, for any e, e' <E {±}< A self-adjoint operator Q acting on is 
said to be P° -trace class Q2] if Q £ 6 2 (£) and Q++, Q~ £ 6i(£). We then 
define its P° -trace as 

trpo(Q) =tr(Q— )+tr(Q++). 

pO 

The space of P° -trace class operators on fj will be denoted by & 1 ~ (Sj). We refer 
to [12] where important properties of this generalization of the trace functional 
are provided. 

2.1 Ultraviolet regularization 

It is well-known that in Quantum Electrodynamics a cut-off is mandatory [Slll7j. 
There are two sources of divergence in the Bogoliubov-Dirac-Fock model. The 
first is the negative continuous spectrum of the Dirac operator, which is cured 
by the subtraction of the (infinite) energy of the Dirac sea, as explained above. 
The second source of divergence is the rather slow growth of the Dirac operator 
for large momenta: D° only behaves linearly in p at infinitjU. 

2 Notice a model similar to the reduced-BDF theory was recently studied for non-relativistic 
crystals in the presence of defects [2] , in which case a cut-off is not necessary because of the 
presence of the Laplacian instead of D°. 
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This can be cured by imposing a sharp cut-off on the space, i.e. by replacing 
£ 2 (M 3 ,C 4 ) by its subspace 

£a:={/g£ 2 (M 3 ,C 4 ) I supp(/)CB(0,A)}. (9) 

Notice D°f)A C f)A. This simple approach was chosen in previous works [T2]— 

om 

However, when looking at decay properties of the electronic density it might 
be more adapted to instead increase the growth of the Dirac operator at infinity. 
This means we replace D° by the operator 

2 * 



Df(p):=(a.p + /3)(l + c(^-)) (10) 



where £ : [0, oo) i— > [0, oo) grows fast enough at infinity. The operator D** is 
self-adjoint on Sj — L 2 (M. 3 , C 4 ) with domain 

I 12 \ \ V 2 



V{D C ) := j/ G L 2 (M 3 ,C 4 ) | + \p\ C ) ) f(p) G L 2 (R 3 , C 4 ) j . 
We remark that the case of the sharp cut-off ([9]) formally corresponds to 

c(*)=(° if t jr'- l; (id 

1 otherwise. 



-oo 



In this work, we shall consider both cases (O and (|10|). We assume through- 
out the whole paper that 

• either f) = Sj\ and £ = (or equivalently C given by (fTTj) ); 

• or ,f) = L 2 (R 3 ,C 4 ) and £ satisfies the following properties: 

C G C 3 ([0, oo)) is non-decreasing and ((0) = 0, (12) 
C(x) > ex e/2 t(x > 1) for some e > 0, (13) 

(l + \x\ p ) C ip) (x) <C(l + ((x)) for p = 1,2, 3. (14) 

Many of our results will be true under weaker assumptions on £ but we shall 
restrict ourselves to (|T2"T) -(fT4" |) for simplicity. We notice that under these as- 
sumptions, the spectrum of is the same as the one of D°: 

ct(D c ) = (-oo;-l] U [l;oo). 

Also the negative spectral projector of is the same as the one of D°: 

P- =X(-oc,o]OD°)=X(-oc,o]OD C )- 

In the whole paper, we shall consider perturbations of of the form D** + 
p* | • | _1 where p belongs to the so-called Coulomb space 

C := {pE S'(R 3 ) | D(p, p) < oo} (15) 
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(16) 



where 

D(f,g)=A-K f \k\- 2 J{k)9(k)dk. 

Jr 3 

Notice the dual space of C is the Beppo-Levi space 

C := {V G L 6 (M 3 ) | W € L 2 (R 3 )} . 

Lemma 1. We assume that fj = f) a and ( = 0, or that fj = L 2 (R 3 , C 4 ) and £ 
satisfies (|12p - (|14p . For any p £ C, the operator + p * \ ■ | _1 defined on the 
same domain as is self-adjoint and satisfies: 



p * 



I — X\ 



= tr eBB (D c ) = (-c»,-l]U[l, 



Proof. We denote V := p 



We have I / |-D'>| 1 is in 6 6 (.?j), hence is 



compact. This is because we can use the Kato-Seiler-Simon inequality (see [24] 
and [551 Thm 4.1]) 



ll/HV) 5 0r)|| e < (2tt)-^ \\g\\ LP( 



\\Lp( 



(17) 



Vp>2, 
and obtain 

Il»(r») IH^Orlz-p.) < C = C \p\ c • (18) 
Lemma [T] is then an application of a criterion by Weyl [221 Sec. XIII. 4]. □ 



M^HIe 6 (*)< C l 



2.2 Definition of the reduced-BDF energy 

We recall that f) = L 2 (R 3 , C 4 ) or f) — Sja depending on the chosen cut-off. We 
need to provide a correct setting for the rBDF energy. When fj = f)A, this was 
done in [H]— [IS]- When fj = L 2 (M 3 ,C 4 ), this is done similarly to the crystal 
case studied in [3]. We introduce the following Banach space: 

Q ■= |Q £ Gift) | Q* = Q, \D C \ 1/2 Q e 6 2 (fj), 

|i)«| 1 /2Q++|i)C|i/2 g 6l(j3); |1?C|1/2q ]i?C|l/3 G 6i(fl)| (19) 



with associated norm 



|^| 1 / 2 Q++|^| 



1/2 



©l(fl) 

|D C | 1/2 Q— |D C | 1/2 



(20) 



We notice that when f) = 55a and ( — 0, one has Q = & 1 (S)a) as chosen 

pO 

in [12] — [15]. In the general case, we only have Q C 6 1 ~(fj). We recall that 
61(f)) is the dual of the space of compact operators acting on f). Hence 61(f)) 
can be endowed with the associated weak-* topology where A n — 1 A in 61(f)) 
means that tr(il„if) — * tr(A?T) for any compact operator K. Together with 
the fact that 62(f)) is a Hilbert space, this defines a weak topology on Q. 
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(21) 



We also introduce the following convex subset of Q: 

K- ■= {Q G Q | -P° < Q < Pi} 

which is the closed convex hull of states of the form Q = P — P°_ S Q where P 
is an orthogonal projector acting on $y. It is clear that K is closed both for the 
strong and the weak-* topology of Q. As we shall see, the reduced BDF energy 
will be coercive and weakly lower semi-continuous on JC. 

Besides, the number tr P o (Q) can be interpreted as the charge of the system 

measured with respect to that of the unperturbed Dirac sea P°, see [H] — |16j . 
Note that the constraint -P° < Q < P£ in flUJ) is indeed equivalent [THT2"] to 
the inequality 

< Q 2 < Q++ - Q— (22) 

and implies in particular that Q ++ > and Q < for any Q £ JC. 

We need to define the density pq of any state Q £ Q. When 9) = Sj\, 
this is easy as any Q e Q has a smooth kernel Q(x,y) (this is because the 
Fourier transform Q{p,q) € L 2 (B(0, A) 2 )). This property was used in [T2"]-[TB] 
to properly define the density of charge. In the case where Sj = L 2 (R 3 , C 4 ) and 
C 7^ 0, this is a bit more involved. The following is similar to |14[ Lemma 1] and 
[U Prop. 1] (we recall that C was defined above in (fT5)) ): 

Proposition 2 (Definition of the density pq for Q £ Q). We assume that 
f) = ft A and C = 0, or that S) = i 2 (R 3 ,C 4 ) and ( satisfies (fl2 ]l -([lg ]> . 

Let Q £ Q. Then QV £ 6f _ (5) for any V £ C . 
constant C (independent of Q and V) such that 

\tx P o{QV)\<C\\Q\\ Q \V\ cl . 

Hence, there exists a continuous linear form Q £ Q<~> 

tx P o(QV)= c>(V,p Q ) c 



Moreover there exists a 



Pq £ C which satisfies 



for any V £ C and any Q £ Q. Eventually when Q £ Q n ©i(^), then 
Pq{x) — tT<£A Q{x,x) where Q{x,y) is the integral kernel of Q. 

The proof of Proposition [5] is given in Section I3~T1 below. 
Let us now define the reduced Bogoliubov-Dirac-Fock (rBDF) energy. In the 
whole paper, we use the notation, for any Q £ Q, 



tr P o (D<Q) := tr (|^| 1 / 2 (Q++ - Q )\D C \ 



1/2 



(23) 



When D^Q £ & 1 (Sj), this coincides with the definition of the generalized trace 
introduced above. The rBDF energy reads: 



£"(Q) = ti'po (D^Q) - aD{ VlPQ ) + -D(p Q ,p Q ) 



(24) 



where we recall that D(-,-) was defined in ((TB|) . In |2^|) . v is an external density 
which will be assumed to belong to i 1 (R 3 ) DC. We use the notation L 3 v = Z. 
The energy is well-defined [12 H3] on the convex set IC. By ([2"2"|) . we have 



tr P o(D^Q)= \D<\ 1 ' 2 Q ++ \D<\ 



1/2 



©i(fl) 



\D^\ X ' 2 Q— \D<\ 



1/2 



(25) 



9 



Together with 

a a 
~uD{v, p Q ) + -D(p Q ,p Q ) > ~-^D{v, v), 

this proves both that is bounded from below on JC, 

VQe/C, e^{Q)>~D{v,v), 

and that it is coercive for the topology of Q. 

Since £" is convex on K. and weakly lower semi-continuous, it has a global 
minimizer Q vac , interpreted as the polarized vacuum in the presence of the 
external field created by the density v. This was remarked in [131 Theorem 3]. 
Assuming that ker(Dg^ c ) — {0} where 

D Q^ = ^ + "(PQvac - ") *l ' 

is the mean field operator, then one can adapt the proof of (T31 Theorem 3] 
to get that Qvac is unique and is a solution of the nonlinear equation Q vac = 
X(-oc,o](DQ mc ) — P_. The charge of the polarized vacuum is —ego where 

?0 = tr p0 (Qvac)- 



When aD{v,v) l l 2 is not too large [T3, Eq. (15)], it was proved that qo = 0. 
However in general electron-positron pairs can appear, giving rise to a charged 
vacuum. When ker(Dg va ) ^ {0}, then £" does not have a unique global mini- 
mizer on K., but it will be proved that qo is anyway a uniquely defined quantity. 



2.3 Existence of minimizers with a charge constraint 

We are interested in the following minimization problem 



E?(q) = inf £" 
QeQ(g) 



where the sector of charge —eq is by definition 

Q(q) := {Q G Q, tr P o (Q) = q} 

and q is any real number. Of course in Physics q € Z but it is convenient to 
allow any real value. It will be proved below that q — > E" (q) is a Lipschitz 
and convex function. Notice that if Q is a global minimizer of £" on Q, then 
qo = tr P a (Q) minimizes q — > E"{q). 

The existence of minimizers to (|26|) is not obvious: although £" is convex 
and weakly lower semi-continuous, and Q(q) is itself a convex set, the linear 
form Q i— ► tr P o (Q) is not weakly continuous. Hence Q(q) is not closed for the 
weak topology. Our main result is the following theorem, whose proof is given 
in Section [321 below. 

Theorem 1 (Existence of atoms and molecules in the reduced BDF model). 
We assume that Sj = Sj\ and ( — 0, or that = L 2 (R 3 ,C 4 ) and £ satisfies 
(THJ) -(Till) - Let bea>Q,ve L 1 (R 3 )nC and denote Z = J R3 v£M. Then there 
exists q rn £ [—00,00) and qM £ [<Zm,°o] such that 
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(*) [?m,9M] is the largest interval on which q — > E"(q) is strictly convex. If 
q M < oo, then E?(q) = E?(q M ) + q- Qm for any q > q M ■ If qM > -oo, then 
E%(q) = E^(q m ) + q m - q for any q < q m ; 

(ii) the interval [q m ,qM\ contains both Z and the unique minimizer qo of q — > 

(Hi) if q £ [q m iQ.u\i then £" has no minimizer in the charge sector Q(q); 

(iv) if q G [q m , 1m], then £" has a minimizer Q in the charge sector Q(q). 
This minimizer is not a priori unique but its associated density pn is uniquely 
determined. It is radially symmetric if v is radially symmetric. The operator Q 
satisfies the self- consistent equation 



Q + P°= X(-oo^) (Dq) + 5, 
D Q =Dt + a(p Q -v)*\-\- 1 , 



(27) 



where fi S [—1,1] is a Lagrange multiplier associated with the charge constraint 
and interpreted as a chemical potential, and S satisfies < S < 1 and Ran(<5) C 
ker(Z?Q — fi). If p G (— 1, 1), then 5 has a finite rank. If p G {— 1, 1}, then 6 is 
trace- class. 

Moreover, pq belongs to L 1 (R 3 ) and satisfies 



7 V- Z 
On — Z = t 

R3 l + aB c A (0) 



(28) 



whe 



B C A (0) = - 



z 4 /3 



I f f 2 - TY ^ = tt logA + 0(l) 



if ft = L 2 (R 3 , C 4 ) and ( ^ 0, and 

T z 1 - z i /3 



B°(0) = i 

n Jo 

if Sj = $j\ and £ = 0. 



2 , A 5 

, dz = — log A 

2 3tt S 9tt 



2 log 2 
3tt 



0(1/A 2 ) 



The constant B^(0) is the value at zero of some real function B^ which will 



be defined later, see ([48]) and 

Equation (|28p has an important physical interpretation. Consider for in- 
stance a nucleus of charge eZ in the vacuum, and assume that Z and its 
distribution of charge v are chosen to ensure that there is no pair creation 
from the vacuum, trpo (Q) = q = 0. A sufficient condition is for instance 

an 1 / e 2 11 / 6 D(v, v) 1 / 2 < 1, see [T3] and Lemma ITU By 025), the electrostatic 
potential which will be observed very far away from the nucleus is a p hy S Z/M 
where 

aphys = TT^Iw ' (29) 

This leads to a new definition of the physical coupling constant called charge 
renormalization (recall that a = e 2 ). The above value of the physical charge 
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Figure 1: Schematic representation of the result. 



(HU) is very well-known in QED, see e.g. [H Eq. (8)] and [H Eq. (7.18)]. This 
was already used and interpreted in [13 , in particular in connection with the 
large cut-off limit A — > oo, in the case = 

The renormalized charge is only observed far away from the nucleus. Close to 
it, one will observe a different behavior like the oscillations of the polarization of 
the vacuum pq. See 16J for an interpretation in terms of the Uehling potential. 

Equation (|28[) implies that a minimizer Q in the charge sector q ^ Z is never 
trace-class, as this would imply trpo Q = L 3 pq and contradict (|28[) . This shows 
that the generalization of the reduced BDF energy E" to the Banach space Q 
is mandatory, as no minimizer exists in the trace class. The mathematical 
difficulty that a minimizer is not trace-class is well interpreted physically in 
terms of charge renormalization. 

When q = Z, it is in principle possible that a minimizer Q for E v {q) is 
trace-class. We shall not investigate this question in this article. 

2.4 Ionization: an estimate on q m and qu 

In the previous section, we have proved the existence of an interval [q m , 1m] for 
q in which minimizers always exist. We now want to provide an estimate on q m 
and qM- We do that with a specific choice for the cut-off function namely 
C(t) = t, which obviously satisfies our assumptions (fT2|) -(fT4 jl . We give this result 
as an illustration: we believe that a same kind of estimate can be derived for 
other cut-offs. The advantage of this choice is that = (—ia • V + /3) (l — ^) 
is local. Notice in this particular case 

*o>-£*a-± + H£H + (*£). 

Theorem 2 (Estimates on q m and qM when Z > 0). We assume that io = 
L 2 (R 3 ,C 4 ) and = t. There exists universal constants < 8q < 1, > 
and C > such that the following holds. For any < a < olq, for any radial 
function v > in L 1 (R 3 ) n C such that Z = J v > and aD(v, v) < da < 1 and 
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any cut-off A > 4 such that a log A < 1/C, the following estimate holds true: 

(30) 



^ ZalogA+l/A + a£>(y, v) . . n 
t> : ~ — ; : S °m S U = qo, 



Z <q M < 



1-Ca log A 
2Z + C(ZalogA + l/A + a% i/)) 



(31) 

oo such that 



1-Ca log A 

In a nonrelativistic limit in which one takes a — ► 0, A 
a log A — > and z/ fixed, one obtains the usual estimate of [20] 

= g m = Qa < Z < q M < 2Z. 

The proof of Theorem is given in Section [331 An estimate more precise than 
(l3"0f and (j3"Tj) is contained in our proof but we do not state it here. 



3 Proofs 

3.1 Proof of Proposition [2] 

When ft = Sja and ( = 0, Proposition's contained in [TU Lemma 1]. Hence we 
only treat the case Sj = L 2 (R 3 , C 4 ). Consider some Q e Q and V S C'ni°°(lR 3 ). 
We have (QU) ++ = Q ++ yP? + Q+-[KP°]P° and (QV) — = Q~7P° + 
Q h [V,P°]P°. We first give an estimate on the commutator [V,P°]. 

Lemma 3. We assume that = S)a and £ = or S) = L 2 (R 3 ,C 4 ) and £ 
satisfies (fl2|) - (fUj) . We have for all r > 1/2 and allp>2 

W, |||^|-[y,pO]|| epW <C||Vy|| iP(K3) 

where the constant C is independent of Q (hence of A J if t > 1/2 or p> 2. 
Proof of Lemma\^ Using Cauchy's formula, wc infer 



[V,P° 



1 

2^ 
1 

2^ 



D° + in 
1 



i 

"2n 



1 1 

a • (W) 



<7r/ 



"7 



rfr/ 



(recall P° does not depend on £). Hence, by means of 



D° + in 



< 



we obtain 

\\\D<\- T [V,P°_} 



1 



\DC\r(D° +irf) 



a ■ (W) 



dt] 
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Next we use the Kato-Seiler-Simon inequality (|17[) and obtain by |T 



\m- T [V,P°}\\ epm <C\\VV\\ LPm x 



dt] 



Lp(R 3 ) 



y/T+rf 



(32) 



(i + |.|r(i + e |.|7A £ )Vi + l-l 2 + ^ 2 

which allows to conclude. □ 
We consider first for (QV) ++ and use Lemma [3] with p — 2 and r = 1/2, 



|Q+-[V,P°]P 



OidO|| <• r 

< C 



Q + -\D^/ 2 
Q\D C \ 1/2 



e 2 (J5) 



©a(.fi) 



|W| 



£ 2 ( 



Similarly we have 



|Q ++ ^|| ei( ,)< 



Q ++ \D C \ 



1/2 



Si($5) 



6c(fl) 



On the other hand, we have by the Kato-Seiler-Simon inequality (|17p 



< 



i3(ffi3) <C||F|| i6 



where we have used Assumption (JTSJ) on £ and |||ZK(-)| 1 l^sfgs) < 00 ■ Hence 



eoc(A) 

by the critical Sobolev embedding i/ 1 (M 3 



<c||VF| i2R3) = c|y| c/ 

L 6 (R 3 ). As a conclusion, 



|tr(QV0 ++ | < ||(Q^) ++ ||e l(fl ) <C\\V\\ C , \\Q\\ Q 



□ 



The proof is the same for (QV) 
3.2 Proof of Theorem QQ 

Step 1: Existence of a minimizer if some HVZ conditions hold. 

Let us start with the analogue of [TU Lemma 3]. 

Lemma 4. We assume that f) = F>\ and £ = 0, or that Sj = L 2 (R 3 ,C 4 ) and 
C satisfies (fT2" j) -fi ty . Let be a > 0, A > and v € L 1 (M 3 ) fl C. We have the 
following estimate 

\q\-^D(u,u)<E^q)<\q\. (33) 
In particular we get for v — and for any ggR, 

E?(q) = \q\. 

Proof. It suffices to follow the proof of [HI Lemma 3] . □ 
Next we state a result analogous to Q31 Theorem 3] . 
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Theorem 3 (A dissociation criterion). We assume that $} — Jo a and ( = 0, 
or that $j = L 2 (M 3 ,C 4 ) and ( satisfies (fT2)) - (fli)) . Let be a > 0, A > and 
v G L 1 (M 3 ) H C. The following two conditions are equivalent 



(H2) each minimizing sequence (Q n )n>i for E^(q) is precompact in Q and 
converges, up to a subsequence, to a minimizer Q of E"(q). 

When it exists, such a minimizer Q satisfies the self- consistent equation 



where p G [— 1, 1] is a Lagrange multiplier associated with the charge constraint 
and interpreted as a chemical potential, and d is a self-adjoint operator satisfying 
< 5 < 1 andRan((5) G ker(Z?Q — p). The operator S is finite rank if p G (—1, 1) 
and trace-class if p G {— 1, 1}. 

Remark 1. Like in [14j Prop. 8], it can be proved that 



In particular this implies that q 1— > E"(q) is Lipschitz. 

Proof. The proof of Theorem is an adaptation of previous works and it will 
not be detailed here. In the case of the sharp cut-off Sj = Sj\ and ( = 0, 
this is contained in the proof of [IH Theorem 3]. In the smooth cut-off case 
Sj = L 2 (R 3 ,C 4 ) with £ 7^ 0, it suffices to follow the proof given in the crystal 
case in [3]. Notice many commutator estimates proved in It (like i], Lemma 
11]) are derived using the regularity of £ and the fact that its derivatives grow 
at most algebraically as expressed by our assumptions (fT^|) -(fLlll. 

The proof that a minimizer Q satisfies Equation (|34ll is the same as in |13[ 
Theorem 3] and O Proposition 2] . Finally, 5 is finite-rank if fj, < 1 because the 
essential spectrum of Dq is the same as that of D° by Lemma [T] If fj, = 1, let us 
recall [HEI] that Q vac := X(-oo,o)Pq) -Pie 6 2 (£) (see Lemma©. By [H 

p" _ p° 
Lemma 2] , we have Qvac G ©1 - (•£))• Hence we deduce Q — Q vac G 6-^(5]) which 

tells us that Q — Q vac and S are trace-class because they are nonnegative. □ 

Proposition 5. Minimizers of ' E"(q) are not necessarily unique, but the density 
Pq is itself uniquely defined. If v is radially symmetric, then so does pq . 

Proof. Note Q G Q — > £^(Q) is convex but not strictly convex. The term / — > 
D{f, f) is strictly convex but the map Q — > pq is not one-to-one. This, however, 
implies that the density pq of a minimizer is uniquely determined, meaning that 
if Qi and Qi are two minimizers of E"(q), then necessarily pq 1 — pq 2 . 

Next we recall that any unitary matrix U G SU2 can be written U = e~ tSn ' a ' 
where G [0, 2ir) and n is a unit vector in M 3 . There is an onto morphism which 
to any such U associates the rotation Rg^ n in R 3 of angle 6 around the axis n. 
The group SU2 acts on 4-spinors in L 2 (R 3 , C 4 ) as follows: 



(Hi) E?(q) < E?(q') + \q- q'\ for any q' ? q; 



Q + P°= X(-oo,„) (D Q ) + 5, 
D Q =Dt + a{p Q -v)*\.\- 1 , 



(34) 



E»(q)<E»(q') + \q-q'\. 



(35) 
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It is well-known [STJ that the Dirac operator D° is invariant under this action. 
As is equal to D° multiplied by a radial function in the Fourier domain, is 
also invariant. When v is a radial function, we hence have £"{Q) = (U QU -1 ) 
and tr P o (Q) = tr P o (U QU^ 1 ) for any Q £ K, and any U £ SU^- This means that 
if Q is a minimizer for E"(q), then U QU~ X is also a minimizer. As Puqu- 1 ( x ) — 
PQ{Rg l n x) , we deduce by uniqueness that pq is a radial function. □ 

Step 2: The density of a solution is in L 1 . 

We prove the important 

Theorem 4 (The density of a solution is in L 1 ). We assume that fj = Sja 
and C = 0, or that Sj = L 2 (M 3 ,C 4 ) and ( satisfies HUH • £ei 6e a > 0, 
A > 0, i/ e L^M 3 ) n C and denote Z = J R3 v £ M. If Q £ Q(q) satisfies the 
self- consistent equation (|?7|) , i/ierc pq 6 i 1 (R 3 ) arte? 

/ po — Z = q ~ Z r . (36) 

Proof. We shall do more than proving that pq £ L 1 (IR 3 ). Namely, we shall 
provide a precise estimate on |/OQ vac H^ip^) needed for the proof of Theorem [21 
Let Q £ Q{q) satisfying the self-consistent equation Q = X(—oo,/j,) {Dq) ~ 
P°_ + 5 where S is a trace-class self-adjoint operator with Ran(5) C ker(£>Q — p), 
and Dq is the mean-field operator: 

D Q =DC + a{p Q -u)*\-\- 1 . 

Recall that by Lemma[TJ o- css (Dq) — cr 0SS (ZX) = (— oo, — 1] U [l,oo), i.e. that 
o-(Dq)D(—1, 1) contains eigenvalues of finite multiplicity, possibly accumulating 
at —1 or 1. For the sake of simplicity, we shall assume that ^ <t(-Dq). The 
following proof can be adapted if £ cf(Dq) by integrating on a line e + ir\ 
instead of ir/ in the integrals below. We introduce the notation 

Qvac := X(-oo,0](Dq) ~ P°, 7 = Q-Qvac- 

Notice 7 £ 6i(fi). We recall that Q++ = P£Q V acP£, Qvac = P-Q vac P° € 
©i(-$3a)- Hence we have to prove that Pq+- + q-+ belongs to L 1 (R 3 ), which we 
will do by a bootstrap argument on the self-consistent equation. 
We can use Cauchy's formula as in [T2] 

1 f°° ( 1 1 \ 3 

with 

Qk = (-l) fc+1 ^ j^— (v'qJ^-) drj, (38) 
2n J_ 00 D<> +ir) \ w + ir] J 

and where we have used the notation ip'q = (pq — v) * \ ■ \ ~ 1 . By Furry 's Theorem, 
it is known that pq 2 = 0, see [HI page 547]. 
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Lemma 6. Let be < r < 1/2. There exists a universal constant C such that 
the following hold: 

lli£ c rQi|| e2(fl) <cii PQ -Hi c , 



|# c | 1/2+r Q; 



li^rai^ri 



e 3 /2(S) 



<c|a»o-Hc> \W dC \Q3 



ls 6/6 (fl) 



< c Wpq - v \\c + a Wpq - u Wc + ' 



< C\p Q - v\ c , 

; Wpq ~ v tc \ 

dist(<7(D Q ),0) i 



Proof. By the residuum formula, we have = Q 1 = 0. On the other hand, 



2tt J^ 00 DC+ ir) 
Hence using 



DC + ir] 



pU 1 /-0O pU 

|£)C|i/a+r 



p o 



D(r,) • 



D(r7)V2- 



where D(r?) := yl + 7 7 2 i an d using also Lemma [31 we obtain 

di] 



C|pQ 



since ||V(/?q| 



Ipq 



We then turn to Q2, inserting 1 = P° + P 1 ^ m fl38|l . We first notice that by 
the residuum formula, 



roo pO / P° \ /" 



pO 



pO 



^ DC + irj V^ Q DC + 277 



d?7 = 0. 



For the other terms, we write for instance 



,00 pO ^ pO ^ pO 

7.^ DC + ir) lfiQ DC + Z77 Vq DC + ?>; ^ 

/oo pO pO pO 

^^^'^DTt^^DCT^^ (39) 

as we did before. We recall that (p'g S L 6 (E 3 ) by the Sobolev inequality. Hence, 
by PU) 



^DT 



177 



< 



6 6 ($5) 



(40) 



Using again Lemma [3j we obtain 



DC 



it] 



pO pO 

+ ir\ ft* DC + in 



< 



< 



C\ P Q 



/2-3/4 



e 2 {S}) 
drj 



, 1 

^ Q D~C~+ir) 



drj 



p (r) )3/4- T /2 
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The proof is the same for all the other terms. 

The same method can be applied to Q3. Let us treat for instance 



A:=\D<\ / 

J — ( 



PS. 



ZX + irj^Q'Di + it] y 4 ZX + ir\ 
00 P° + \D<\ r , DOl P°- 



P° 



dr) 



Applying the above method with 



1 

D^ + ii] \ ^ Q ~D< + ir] 



dr). (41) 



[<p> Q ,P°_] \D<\-^ 
by Lemma |3l we obtain 



< 



6 2 (f>) 



clival 



L 2 (R3) 



\\ 66M V<C\\p Q -v\\ c J^ (1+??2)5/8 



dr) 



(42) 



The argument is of course the same for all the other terms. 

Finally, we expand further Q' A to the 6th order: Q' 4 = Q4 + Q$ + Q' 6 where 
Q4 and Q5 are given by (|38|) and 



Q'e 



1 

'2?r 



i D^ + iri V>Q J D Q + irj \J Q D< +ir] J 
On the one hand, we know that \Dq + irj\ > dist(cr(Z?Q), 0)), and therefore, 

\\{D Q + <&st{a(D Q ),Qy\ 
On the other hand, we can use (FIT)]) and 



fQ- 



dn. 



ZX + ir\ ' 



e 6 (Sj) 



< c \\\D (-)r II £6^3) Ik'olLe^) 



to estimate |ZX| T Qg|ZX|' r . The terms Q4 and Q5 are treated like Q2 and Q3. □ 
Lemma 7. Let be < t < 1/2. There exists a universal constant C such that 



(43) 
(44) 



IW^Qf^rW^ < C (\p Q - vf c + a 2 \\p Q v\ 

Ill^rO^l^rlle,^ < C (\\p Q - uf c + a 2 \\ PQ v\ 
Proof. Consider the operator 

D{t) :=fl c +t kiM^ , 

Wpq ~ v \c 

Since pq — v G C, we can use (|18l) to deduce that there exists a universal constant 
<o > such that \D(t)\ > 1/2 for all t e [— £ , to]- Next we introduce 

Q(t) ■■= X^-MiD®) - P°. 
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We can write as before 

= E „ tL ^ Qfe + „ * 4 (*) 

k =i \\PQ ~ flic \PQ ~ v \c 
where Qk are denned as above and with this time 

Following the method of Lemma [6] and using |-D(£)| > 1/2, we can prove that 
\\\D c \ T Q , MD C \ T \\ eaW <C[pQ-u\*, 

WlD^Q'MD^W^ < C (lp Q - uf c + a 2 \\p Q - uf c ) . 
Next, the estimates of Lemma [5] imply that 

for all t € [— to, to]. But as Q(t) is a difference of two projectors, we have 
Q{tf = Q(t)++ - Q{t)— £ Si(5). Thus 

|||i> c l T 0(*) ++ |£> c l T || ei(J5) + \\\D^Q(t)-\D^\\ 6im < a 

Finally 

— ^—^Qt + + — ^—^Qt + = QW ++ i —^Q' i {t) ++ 

\\PQ ~ v\c WQ - v \c iPQ - Hie 

which gives the result when applied to t — to and —to. □ 
Lemma 8. Let be 2 < p < 6. There exists a universal constant C such that 

\\P(Q 3 )+- * I ' TIl^) + \\P(Q»)~+ * I ' I _1 |Il P (r3) <C\\pQ-vf c . 

Proof. By LemmaH Q 3 \D<\ £ 6 6/5 (ij), hence Q 3 \D^\ e & q (Sj) for all q > 6/5 
and 

\\Qt-\D C \\\ 6qm < |0a|I? C |||e fW <C\\p Q -vf c . (45) 
Let us choose a test function V in the Schwartz class. We have 

|tr((Q 3 ) + ^)| = |tr((Q 3 )+-P^P«)| = | tr ((Q 3 )+- \D<\ \D<\- l [P°,V]) \ 

<lW3) + -|iJ C |||e f(J9) ||l^|- 1 [i*,V]|| eff(J9) (46) 

for all q > 6/5 and q' = q/(q— 1). Then we use Lemma [3] which tells us that 

\\\D^[Ptv}\\ eA ^<c\\vv\\ Lt , m 

provided q' > 2. Finally by the Sobolev inequality and Riesz operator theory 
|W||^ (R3) < C\\D 2 V\\ lp , (r3) < C'\\AV\\ LP * {R3) 
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for p* = 3q'/(3 + q'), 2 < q' < 6. Summarizing, by (ji5)) and igBJI , 
|tr((Q 3 ) + -^)| <C\\p Q -vf c \\AV\\ LP , {R3) 



for any 6/5 < p* < 2. By duality, this proves that for any 2 < p < 6 

\\p(q 3 )+- * i • r 1 IL P (R3) < c ipq - hi c ■ 



□ 

Lemma 9. Let be 3 < p < oo. There exists a universal constant C such that 
lhQ3)±± * I • T 1 1| L p (R 3) < C (|1pq - vf c + a 2 \p Q - u[ty , 

-111 ^r<(\\„ ,.l|4 , „.«„ ,.i|5 , „2 IPO — Hie I 



PQi * I • I L,™ <c Ipq - Hie + 01 \PQ ~ HI 



| iP(K3) ^ - ^IFW -HC < "IPW -IIC • - dist( ( rp Q ),0) / 
Proof. We argue as above, taking some V in the Schwartz class. We have 

itr(Q++v)i < m^rg+ + i^rii ei(fl) n^r^i^riic^) 

for any q > 1 and r < 1/2. Then by the Kato-Seiler-Simon inequality (fTT|) 

li^rvi^i-|| e€(a) < c||£(.)- 2 i L9(R3) mi i9(R3) 

which makes sense as soon as q > 3 and 1/2 — t is small enough. The rest 
follows from the Sobolev embedding like in the proof of Lemma [HI □ 

We now estimate PQ vac using the self-consistent equation. First we recall 
that pQ 2 = and that Qi = (Qi) + ~ + (Qi) h can be explicitly computed [T2] 
yielding 

p^Z{k) = -aB{(k) {p^Z{k)+p^{k)-v{k)) 

+ « 4 P^(fc) + « 3 P^p (*) + a 3 p^Tl (fc) + a 3 /5^p (fc) + (fc) (47) 

with 

1 f (£ + k/2)-(£-k/2) + l-E(l + k/2)E(£-k/2) 
E(£ + k/2)E(£-k/2) X 

' ^rd£ (48) 



£(l + fc/2) (l + C (^^) fc/2) (l + C (— 

when = L 2 (R 3 ,C 4 ) and C satisfies (PD-flHI), and 

1 / (£ + fc/2) ■ (I - fc/2) + !-£(* + fc/2)g(l - fc/2) 

Al J 7r 2 |fc| 2 J\i+k/2\<A, E(e + k/2)E(£- k/2)(E(£ + fc/2) + fc/2)) 

(49) 

when F) — Sja and £ = 0. Notice that in both cases B^ is a radial function. 
Also B A has its support in B(0, 2A). 



20 



Remark 2. There is a small mistake in the domain of integration of the defi- 
nition of B A in O Eq. (40)]. This does not change the analysis of [T^] but is 
important for the present study. 

Many properties of B A and B A are given in Appendix. Let us define b A by 



b{(k) = (2nf/' 



\ + aB A (k) 

2*3 



(50) 



In both cases f) = ¥)a with £ = 0, and S) — L (R, C ) with £ satisfying |T 
(fT4|) . we prove in Appendix that b A is a smooth function belonging to L 1 



see Propositions [17] and [18] In the rest of the proof, we use the notation 



Ia 



R3 



\b A (x)\dx < 



Equation (|47[) can be rewritten as 



(51) 



PQv 



(V - p 7 - pi) + pi + p 2 - b\ * p 2 , 



Pi = a4 PQi + « 3 Pq 3 ++ + G L^R 3 ), p 2 = a 3 

By Lemma [3 Lemma [9] and ([51"]) 

(p 2 - & A * p 2 ) * | " 



(52) 



a Pq;+- 



L 4 (K 3 



<C(1+/ A ) a 3 ||pq-HIc 



(pi - &x*pi) * i • r 



L 4 (R3) 



< C(l + /a) a 3 Ipq - Hie + « 5 Ipq - vf c + 



\\PQ 



dist(<7(D Q ),0) 



so that 



Il 4 (r3) 



<c(i+/ A )^||(^-p 7 )H-r 1 | L 4 (R 3 ) + 

i 3 ii ii 3 , s ii ii 5 , 6 Ipq ~ v \c 

+ a\\ PQ -v\\ c +a In-^+a—p— 
As p 7 , i/ e C n L^R 3 ), we have 



(53) 



l^-Ar)*!-! -1 ! 



L 4 ( 



< oo 



but we do not provide a precise estimate at this point. Now we can use the 
information that <pg € L 4 (R 3 ) to estimate (Q3)" 1 and (Q3) h , using 

C „ , „ 



For any fixed < r < 1/2, this gives an estimate of the form 

\\\D^Q+-\D<r\\ 6i(fj) < C lp Q - V \ c y Q \\l HR3) . (54) 
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Inserting in (f52")) . we are led to 



+ a C(l + 7a) ||pq - Hie || ( ^Q|| i 4 (R 3) 

(55) 

where C is independent of A. As a conclusion, pQ vac hence /3q belong to L 1 (IR 3 ). 

Let us turn to the proof of We deduce from the previous analysis that 
PQi e i 1 (K 3 ) (whereas in general Qi 6i(£)a)) and that 

J PQ = J PQ+++Q— + a J Pq+-+Q^+ =q + a J PQt'+QT + 

since we know that 7, Q3, Q§ and Q$ are all trace-class and that J P{k)+- — 
tr(P±K P®) = for any trace-class operator K. Now 

PqM = -B A (0)(^(0) - u(0)) = -Si(0)(p5(0) - P(0)) 
which leads to 



^ ( °? (q-Z) and /( PQ - y ) = g ~f . 
l + aS A (0) J W l + a£ A (0) 



This ends the proof of Theorem HJ □ 

Corollary 1. Let Q be a minimizer for E"(q) as in Theorem^ If q < Z (resp. 
q > Z) then o(Dq) contains an infinite sequence of eigenvalues converging to 1 
(resp. to —1). 

Proof. This is a simple adaptation of the proof of [2J Thm A. 12]. □ 

Step 3: Properties of q — > E"(q) and definition of q m and qM- 

As Q G Q — > £^ (Q) is convex, the map q — > £^ (g) is also convex. We then 
define I = {q G M | (HI) holds}, where (HI) is defined in Theorem[3] Thus, for 
any q G I, there exists a Q G Q(q) such that £ r "(Q) = E"(q), by TheoremE] We 
introduce the following convex real functions f~(q) ■= E"(q) — q and f + (q) ■= 
E"(q) + q. By (|35|) and (|33|) . / _ is nonincreasing and bounded from below, 
f + is nondecreasing and bounded from below. Notice lim^oo / + (q) — 00 and 
limq^-oo f~(q) = 00. Define now qM such that f~ is decreasing on (— 00, qM) 
and constant on [<7m,oo) (let qM = 00 if /~ is decreasing), and q m such that 
f + is increasing on (g m ,oo) and constant on [— 00, q m ) (let q m = —00 if / + is 
increasing). Remark q rn < qM- Next we have 



qel 



Vq'>q, E»(q)<E»{q')+q' -q 
Vq'<q, E?(q)<E?(q')+q-q' 

Vq'>q, f+(q)< f+(q>) 
Vq><q, f-(q)<f-(q') 

q G fern, 00) 
q G (-00, q M ] 



and therefore 7 = [q m , qM] 
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Step 4: The interval [q m ,QM] contains both qo and Z. 

Assume now that q$ satisfies E^(go) — min ge u_Ey (g). Then 

E?(q )<E»(q')<E»(q') + \q -q'\ 

for any q' ^ go, ie. qo satisfies (HI). Hence qo 6 / = [g TO , gjvf]. 

Let us now prove that Z = J" i/ also belongs to I — [q m ,QM]- We use 
classical ideas already used for the reduced Hartree-Fock theory [26] . Assume 
first Z > qM- Since G /, there exists a minimizer Qm in the charge sector 
Q{qu)- By Theorem^ Qm satisfies the self-consistent equation 

Qm+P- =X(-oo, M ) (-°Qm)+^ 

for some /j, € [—1,1]. By Corollary [TJ ct(Dq m ) contains an infinite sequence 
of eigenvalues converging to 1. Since tr P o [ X (-oo,o) (Dq m ) — P_] is known to 
be finite and S is finite rank, we deduce that /i < 1. Hence there exists an 
eigenvalue A € (/z, 1) of Dq m with eigenfunction % € $)\ which is not filled. 
Notice Qm + t\x){x\ £ Q(qM + 1) for i e [0, 1]. Let us then compute, 

E»(q M +t)< £»{Q M + t\ X )(x\) = E»(q M ) + t(D QuX ,x) + ^~D(\x\ 2 , \x\ 2 ) 
or equivalently 

f-(q M + t)< f-(q M ) + t(X - 1) + 0(t 2 ) 

which contradicts the definition of q^j. 

Assume now Z < q m and consider a minimizer Q m for E"(q m ). By the same 
arguments, it satisfies the self-consistent equation 

Q m + P°_ = X(-oo lP ')(A? m )+* 

for some // > — 1 and the spectrum o~(DQ m ) contains an infinite sequence of 
eigenvalues converging to -1. Thus there is an eigenvalue A' G (— which is 
completely filled, with eigenfunction x' € Computing £^(Q m — t\x'){x'\) 
and noticing Q m — i|x')(x'| <E Q(q m — t) for any t € [0, 1], we obtain 

f + (q m -t)< f + (q m ) - t(A' + 1) + 0(t 2 ) 
which contradicts the definition of q m . 

Step 5: Characterization of [q m ,qM]- 

Lemma 10. Assume that q\ ^ q2 are such that both E"(qi) and E"{q%) admit 
a minimizer. Then 

Vie (0,1), E^{tq 1 + {\-t)q 2 )<tE^(q 1 ) + {\-t)E^{q 2 ). (56) 

As a consequence, 

I- [9m, 9m] *s fie largest interval on which q — > E"(q) is strictly convex; 
2. 9o = argmin^i?^ is uniquely defined; 
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3. no minimizer exists for E"(q) when q is outside [q m ,QM]j 

Proof. Assume that Q\ and Q 2 are two minimizcrs of respectively E"(qi) and 
E?(q 2 ), with qi ^ q 2 . Then by ([28]), / p Ql ^ J Pq 2 , hence p Ql ^ p Q2 . Hence, 
for any t € (0, 1), 

E?(t qi + (1 - t)q 2 ) < £»{tQ 1 + (1 - t)Q 2 ) < t £ r "(Qi) + (1 - tJ^COa) 

= t££(9i) + (l-t)E»(q 2 ) 

where we have used the strict convexity of / 1— > £>(/, /). 

Inequality (j5l?)) shows that q — > E"(q) is strictly convex on J = [g m ,QM], 
since minimizers are known to exist for any q £ I. But 9 — > E"(q) is linear 
outside I and therefore J is the largest interval on which q — > E"(q) is strictly 
convex. The global minimizer go of on K thus on / is unique. 

Eventually, we prove that no minimizer exist for E"(q) when q £ [q m ,qM\- 
If q > qu provides a minimizer, then since a minimizer exists for E"(c[m), (156p 
applied for qu and g contradicts the fact that E"(-) is linear on [qM, 00). □ 

3.3 Proof of Theorem H 

If we assume £(t) = t, the function can be studied more carefully as explained 
in Appendix. In this case, one can prove that 

< aB{(0) < 2/(37r)alogA 
Li(K3) - i- a B c A (0) ~ 1 - 2/(37r)alogA' 

when A > 4 and 2/(37r)alogA < 1, see Proposition [5T] For the sake of simplic- 
ity, we shall use the following notation in the whole proof 

9 := a-K l ' 6 2 xl l 6 D{v, v) 1 ' 2 

and we will assume that 9 < 1. Later on we shall also assume that a, I and 9 
are small enough but we postpone this to the end of the proof and rather give 
precise estimates before. 

Step 1: A priori estimates. 

Lemma 11. Assume that Q G Q(q) is a minimizer for E"(q), for some q £ 
[IttdIm]- Then we have 

\\PQ-v\\ c <\Mc- ( 57 ) 

If moreover 9 := om l l % 2 lx l % D(v, v) 1 / 2 < 1, then \Dq\ > 1-0, hence g cr(D Q ) 
and, denoting Q vac = X{-oo,o){Dq) — P_ we have 

tr P o (Q vac ) = 0. (58) 
Proof of Lemma [771 We have by (|3"3"|) 

tr P o (D C Q) + ^D{p Q -u,p Q -v)<^D{v,v) + \q\. 
Introducing q + = tr P o (Q ++ ) > and q~ = — tr P o (Q ) > 0, we have 
tv P o_(D<Q) = tr(|£K| 1 /2(Q++ _ Q—)\D C \ 1/2 ) > q+ + q~ > \q\, 



I A = 
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hence (|5T|) follows. Following [T2J p. 4495], we have the operator inequality 



(PQ -V)* — 



<k\\v\\ c \D°\ <kH c |£»C| 



with k = 7r 1 / 6 2 11 / 6 . Hence 



|A,| >(!-<*« \v\\ c )\Dt\>l 



(59) 



The proof that tr P o (Q vac ) = is the same as in [12 [T3]: considering P(t) = 

X(-oo,o)(D c +at(pQ-u*\-\~ 1 )), we have by [H Lemma 2] that tr P o (P(t)-P°) = 

tr(P(t) — P") 3 is for all t S [0, 1] an integer which varies continuously with 
respect to t, hence, it is equal to for all t G [0, 1]. □ 

For the rest of the proof, we work under the assumptions of Theorem 
namely we assume that v is a radial and positive function in L 1 (R 3 ) n C such 
that an 1 / 6 2 n / 6 D(v, v) 1 ! 2 < 1 and Z = J v > 0. Let Q be a minimizer for 
E"(q), q £ [qm,qAi]- It solves the self-consistent equation 



Q ■= Qvac + 7, 



Qv 



X{-ocM D Q)- P -- 



(60) 



By Lemma HT1 7 is either > if q > or < if q < 0. It satisfies IMIg/^) = \q\- 
As pq is radial by Proposition the operator Dq is invariant under the action 
of SU2 introduced in the proof of Proposition [5] In particular, we deduce that 
UQvucU^ 1 = Qvac for any U G SU2- Hence PQ mc is also a radial function. 
Therefore p 1 = pq — PQ mc is radially symmetric. 

Lemma 12. Assume that Q G Q(q) is a minimizer for E" (q) for some q S 
[imjQM]) decomposed as in ((6*U|) . and that 9 :— onr 1 / 6 2 11 / 6 D(i>,v) 1 / 2 < 1. Let 
be < t < 1/2. There exists a constant C > (depending only on t) such that 



||l^ c l r Qv 



< 



ce 



IMQvacI < 



c 



1 - log(l - 9) 



1 - log(l - 9) 1 



\PQ 



(61) 



(62) 



Proof. We have 
Hence 

||l^ c l T Qv 



2vr 



D 1 ^ + ir]' fQ D Q + in 



—drj. 



Is 6 (fl) 



< 

< C0 

< C0 



v/(i-^) 2 + ^ 2 

drj 



\D Q {P)\ T , 

(\d<( p )\ 2 + y?yi^ Q 

\D C (P)\ T 



- 9) 2 + ?f 



(|^C(p)| 2 +7 ? 2 )l/2 



6 6 (.s) 



< 



C0 



/-00 S(7 ? )l/2-r^( 1 _0)2 +?? 2 " 1 _ log(l - 9) 

by ([5T|) . (|5"9")l and the Kato-Seiler-Simon inequality ([T7|) . 
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Notice |x|<5vac = jfy • xQvac and ar|a;| 1 is bounded on $). Hence for (f6"2"| it 
suffices to prove that XfeQ vac is a bounded operator for any k = f ..3. We write 



x k (. 



Notice 



2tt 



1 



1 



Xk, 



DC + it] 



Xk, 



1 



DC + irj 



1 



f 



f 



D Q +ir) DC +irj Xk<fQ D 



Q+iV 



drj. 



DC + ir] 



D c + irj] 



1 



1 



DC + ir] DC + irj 



(d Pk DC) 



DC 



If] 



Clearly d Pk DC = a k (l + \p\ 2 / A 2 ) + 2p k /A 2 {a -p + 0), hence 



H 2 M 2 



< c. 



and by dSSJ, <0DJ) and (57|) 





f 


J — oo 





1 



< 



< 



ch) 



— HMIc- 



f - log(l - 0) 

Since pq and v are radial, we have by Newton's theorem 

\PQ - v \(v) 



\x k f' Q {x)\ < \x\\<p' Q (x)\ < \x\ 



R3 \x-y\ 



dy < \p Q - v\ 



hence 



dTTwi XWq ~d^~ 



IT] 



-drj 



< 



C 



1 - log(l 



\PQ 



This ends the proof of Lemma [TSJ 
Lemma 13. We have 

hDC 



< 



□ 



(63) 



6i(D) -1-61' 

Proof. Assume for instance that q > and 7 > 0. By the self-consistent equa- 
tion (j6*fj|) . we have 7-Dq > and \.v{^Dq) < tr(7) = 9. Hence IItA^Is^) < q 
We then write 

7^q = 7^ (1 + osgn(D c )|D c |- V«) • 

We now use that 



i S6(fl) <MpQ-Hl c < K Hc<- 



IP c rv Q | 

by [331 p. 4495] and (|57|. so that 1 + asgn(L> c )|L> c |~y Q is invertible and 
(l + aBgn(I*)|lK|-yo)" 



< 



1 - 



This gives the result. 



□ 
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Lemma 14. There exists a universal constant C such that 
(l-C(l + / A ) 3 a0 2 ) / | PQvac | 

< (J A + C(l + i A ) 3 a0 2 )(^ + |g|) + C(l + /a) f^ 2 + I . (64) 



1 - 

Proof. By (|5"5j) . we have 

l|PQ va Jlz,i ( M 3) < W + M) + c (! + j a) + Y?e) 

+ C(l + /A) 3 a 3 ||^||^ 4(R 3 ) IHIc- (65) 

Notice for any p 

\\P*\- \~\*^) ^ G \p\c 2 IpIli 2 (r3) (66) 
which is proved by writing 

\P*\'\ ||l,«(R3) - CjPWIkl 2 |Il4/3 (R 3 ) 

^ |plii(lK3) |||fe| l £ 4/8(B(o,r)) + l^ic l|W IL*(lEl3\B(0,r)) 

and optimizing in r. Using (|66p for p = pq — v and (|57p. we get 

1^11^3) <c He/ IpoI + "<c He (^ + M + / IpqvJ)- 

Inserting this in (jB3|) yields the result. □ 
Step 2: Lieb's argument. 

We now use ideas from Lieb [20] to obtain a bound on qu- We denote by Q 
a minimizer for E IJ (qM) which exists by Theorem Q] As qu > Z > 0, we can 
decompose Q as in (jBTif : 

Q = X(-oo,,o Pq) - P° + 6 = Q vac + 7 (67) 
with 7 > 0. Using that (D Q - 1)7 < due to ([57|. we infer 

> tr(\x\(D Q - 1) 7 ) = tx(\x\(Dt - 1)7) + a [ \x\<p' Q (x) Pl (x)dx. (68) 

JR 3 

Lemma 15. There exists a universal constant C such that 



f n im( i\ \ ^ CqM ( 1 1 aqM + a ^ 
tr(N(D« - 1)7) > -— 9 [j + i _ bg(1 . 



The proof of Lemma [T5] will be given at the end of this section. Now we 
assume that a, I\ and 9 are all small enough. Then (|69|) becomes 

(70) 



tr{\x\(D ( - 1)7) > -Cq M (j + aq M + aZ + d^j 
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and becomes 

\p Qv J < C(J A + «# 2 )(Z + 9m) + C {a6 2 + 4 ) . (71) 

To estimate the second term in (|68p , we write 

on ii | ii 

\xW Q {x)p 1 {x)dx = // -p 7 {x)p 7 {y)dxdy 

JJ R 6 2\x - y\ 

\X\{PQ -UM^X) 

F - V\ 
and notice 

W + \ y \ Pl (x)p,(y)dxdy>4f- (72) 



since p 1 > and |a; — y| < + |y|- Using Newton's theorem we infer 

\x\i/(y)p y (x) 



F - y\ 



-dxdy<Zq M , (73) 



^MM^M dxrf^ < C(I A +a9 2 )q 2 M +Cq M (I A +a9 2 )Z+C qM (a9 2 + 4 ) 

by (|7ip and since both j/, p 7 and pg vac are radial functions. Collecting estimates 
and using that a, I and 9 are small enough, we obtain the following estimate 

C 

(1- Co: log A)g M < 2{l + Ca\ogk)Z + — + C6. (74) 



The proof for q m is the same, using that in this case 7 < and instead of ([73]) 
- // M V } v) ^ x) dxdy>0. 

JJrs \x-y\ 

Proof of Lemma \75l For the second term of (|68|) . we compute 
tr(|x|(L» c - 1)7) = tr(|x|(|L» c | - 1)7) - 2 tr(|x|D c P°7) (75) 
= tr(|s|(|D*| - 1)7) + 2tr(\x\D^Q V3jc j) 

= tr(|.x|(|£ c | - 1)7) + 2tr([|a;|,P c ]Q vac 7) + 2 tr(\x\Q v , cl D c ) 
where we have used that X(-oo,o]{Dq)j = by IpTF]) . One computes 

(\D<(p)\ l)\x\ + \x\(\DHp)\ 1) = (E(p) l)\x\ + \x\(E(p) 1) 

3 3 

+ £ ^r*b*. NO + P*U]^|p + ^ X>»)M + \*\e(p))pu- 

k=l k=l 

Next we use a result of Lieb [20] which says that 

(E(p)-l)\x\ + \x\(E(p)-l)>0. 

We obtain 



3 r 



QDt(p)\-l)\x\ + \x\(\Dt(p)\-l)>^Yl 



A 2 

k=l 
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and 



Notice 



tr(H(|^|-l) 7 )>-^£ 



fe=l 



E(p)pk 



\DC(p)\ 



ItI^I" 1 ! 



E(p)pk = Pk < A 
|ZX(p)| l + |p| 2 /A 2 " 2 



hence, using (j6"3"|) , we obtain 



tr(N(|^|-l) 7 )>- 



C 



(l-fl)A 



<7Af- 



Let us now estimate the last term of the r.h.s. of (|75|) . Using 
obtain the following estimate: 



and 



|tr(|x|Q vac7 D c )| < 



G 



(aq M +aZ + 8) q M . 



(1 - - log(l - 0)) 
Eventually we estimate the second term of the r.h.s. of (ITS)) . We compute 

3 

ru 2 mi j- ' 

A 2 



(76) 



we 



(77) 



[« - ^- I'll + ^[bl 2 , ki] + X) ^Ib/bbl 2 , k|] 



k=l 

k=l 1 1 fe=l 



[bl 2 ,N] = A_ 2ip .^ 

Mp| , = RPfc + 2p k - - 2p ■ x— + 2 iPkP ■ ^ - t|p| R . 

Hence, using Hardy's inequality which tell us that |p| _1 |a;| _1 is a bounded op- 
erator on fj, we easily deduce that 

[\ x \,Dt]=-± i ±. a -p + A 

where A is an operator satisfying |||D^| _1 A|| < C for a universal constant C 
independent of A. Next we write 



> 



2 

A 2 
C 



1 1 



tr([|a ; |,^]Q vac 7) = -— tr(p|^a-pQ vac7 |p|)+tr((^)- 1 AQ vac 7^) 



A 4 /3 

> -C 



\P\ 



^2/3^ vai 



qM 



-C\\Q V 



(l-0)log(l- 



(78) 



by Lemma [TBI and Lemma [T^] with r = and r = 1/3. Inserting (ITB^l . (|77p and 
(|78p in Formula ([75|l . we obtain (|6"9"|) . This ends the proof of Lemma fT51 □ 
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A Study of the function 6^ 

This appendix is devoted to the decay properties of &£, for the different cut-offs 
chosen in this article. The function bj^ plays an important role in the model as 
it can be interpreted as the linear response of the vacuum in the presence of an 
external field, as shown by Formula (|52p. We recall that 

l + aB c A {k) 



b{(k) = (2^/ \ ^ (79) 



A.l Study of 6° when 9) = S)a and ( = 

We start with the sharp cut-off case £ = and Sj = In this case 

B o m = L_ f (g + fc/2) ■ (g - V2) + 1 - g(g + k/2)E(q - fc/2) 

Al J 7r 2 |fc| 2 7k+fe/2|<A, E(q + k/2)E(q-k/2)(E(q + k/2) + E(q-k/2)) 9 

| (j— fc/2|<A 

is defined for |fc| < 2 A. Following [33], for any g £ I 3 we introduce as new 
variables the azimuth angle ip around the axis parallel to k and 

v = (E(q + k/2)-E(q-k/2))/2, 
w = (E(q + k/2) + E(q-k/2))/2. 

Then integrating over {q £ R 3 | \q + k/2\ < A, \q — k/2\ < A} is easily shown 
to be equivalent to integrate over the new variables (u,v,ip) £ K x R x [0, 2tt) 
with the three conditions 



1 < v + w < a/I + A 2 , 1 < w - v < \/l + A 2 , (82) 
VTTPf74< W <v^TA^, M < 1*1 ^5^ 1, (83) 



Eventually (|82[) and ([83]) are equivalent to 



Vl + |A:| 2 /4 < w < Vl + A 2 (84) 



M <^^-i.^- n .ia ^ji^-' j. <k, 

An explicit computation shows that 



1, y/T+A?-w, ^ 



w 2 - |fc| 2 /4- 1 



2 V u> 2 - |/fc| 2 /4 



^V^St? when W <^ A (|fc|) 
Vl + A 2 - w when w > W\(\k\) 



(86) 



where W\[r) := (vT+A 2 + ^/l + (A - r) 2 )/2 is the unique root of the fourth 
order polynomial equation 

r — |fc| lw 2 - |fc| 2 /4- 1 
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in [a/1 + |fc| 2 /4, Vl + A 2 ]. Inserting this in the definition of B\(k) and using 
that dq = (2/\k\)E(q + k/2)E(q - k/2)dvdwdip, see [H Eq. (12)], we find 



n\k\ 3 \ J ^1+^/4 Jo w 

dw / dv l \EJ± . ( 87 ) 

Wa(|As|) Jo w I 



Letting z = \J W w -i^\Jpj± m the first integral and z — 2(\f\ + A 2 — in 
the second, we obtain 

i /•■Ml*|) 7 2_ y iio 



7T J Q (l-Z 2 )(l + |fc| 2 (l-Z 2 )/4) 

\k\ f z ^ z-z 3 /3 



2tt J q - \k\z/2 

where we have defined 



-dz 



W A (r) 2 -r 2 /4-l _ VI + A 2 - y/1 + (A - r) 2 
W A (r) 2 - r 2 /4 ~~ r 



^A(r) 



The first term of (|88| was already present in [22] , whereas the second term was 
ignored by Pauli and Rose. An explicit computation of the integrals in 
yields 

flO(fc) = _L( _ _ 2 (r2 _ 2) ^^ arctanh f ^A(r) 



tt|A;| 3 [ 3 3 V yv Vn/IT^ 
7-3 f 1 + z A(r)\ , 8£(A) 3 f 3r 2 ^ Z A (r) , 44^ 3 



- 2£(A)r 2 - 2 (3£(A) 2 + 1 + r 2 ) Z A (r) + 3£(A)Z A (r) 2 - ^ A (r) 3 j. (89) 

To avoid any further notation, we now see B\ as a function of \k\ . Using Formula 
(f89| , one can prove the 

Proposition 16 (Regularity of £> A ). Let be A > 0. TTie function r i— ► i? A (r) 
extends to a non-negative, C 1 function on R + , which vanishes on [2A, +oo). 
Moreover, it is of class C 3 on [0, 2A]. Eventually, we have 

Bl(0) =Bl = l- ln(A) + 0(1), ^(0) = + ^oj ±) , 



, ,(0) = -—+ O , -— A(o) = — + o — 

dr 2 v ; 157T A^+cc VA 2 / dr 3 v ' 47rA a^+ooVA 2 

WR° d 2 R° A d 3 R° ^A 2 — 1 

S A (2A) = ^(2A) = 0, ^(2A) = — * ^(2A) " 



dr 1 ' ' dr 2 y ' 4tt£;(A) 3 ' dr 3 v ' 87r£(A) £ 
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By Proposition 1 161 B A is a non-negative, continuous function with compact 
support. Therefore, by ((75)) , b A is a smooth function which reads (using the 
inverse Fourier formula for spherically symmetric functions), 



9 p2A 

VxeR 3 \{0}, b° A (x) = - / - 

PI Jo 1 

In particular we get the bound 



+ aB° A (r] 



■ sm(r\x\)rdr. 



K(x)\ < 



16ttA 3 



(90) 



(91) 



which shows that b A S L°°(R 3 ). This becomes after three integrations by parts, 



«W - |L{2aA(BS)"(2A)co 8 (2A W ) + ^3™ 



,-2A 



ar(B° A )^(r) 6a 2 r(B A )' (r)(B° A )" (r) 6a 3 r(B° A )' \rf 



! y(l + aB° A (r)y (l + aB° A (r)f 

Za{B A )"{r) 6a 2 {B° A Y{r) 2 \ 



(1 + aBXW) 2 (l + aB°(r)) 3 y 
which yields by Proposition [TO] 

\b A (x)\ < 



cos(r|a;|)<ir 



(l + aBl{r)Y 
(92) 



(93) 



for some constant C a ^ A depending on a and A. With (f9"Tj) , this proves the 

Proposition 17. Assume io = £)a C — 0- ^e£ &e a > and A > 0. Then 
b° A belongs to L 1 (R 3 ). 

Remark 3. It can be seen that B A (Ar) — * B%^{r) where 

4 r (!^)iog(i)-^# + ^ + 



BUr) 



4 
-ttt 3 



4^9 



when 1 < r < 2, 



■0" 
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2-r\ 



■log I 



when < r < 1. 



The convergence holds in C 2 ([0,2)). Notice the two terms appearing ([55)1 sep- 
arately converge to a function which is not differentiable at r = 1, but there 



is some cancellation occurring. It can also be proved that ||&aIIz,i(r3) * s m deed 
uniformly bounded independently of the cut-off A, but we do not need that in 
this article. 



A. 2 Study of b c A when ft = L 2 { 



and C ^ 



When ft = L 2 (M 3 , C 4 ) and C ^ satisfies fT2") ) -([i"I )) . the same changes of vari- 
ables, followed by t = \/l — u 2 , lead to 



B c A (rk) = tt- 1 



dt 



t(l + \k\ 2 t 2 /4)J 



l-u 1 



l + *(|A|,t,u) 



du, 



(94) 
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where 



1 v 

i$(\k\,t,u) = -(11(11) + v) + r](w — v)) + —Mw + v) — T](w — v)), 

2 Zw 



T)(x) = ( 



x 2 - 1 



\k\u 

v = and w 

2 



kl 1 
4 t 2 ' 



4 ) and that ( ^ satisfies (fT2"j)-(fT 



Proposition 18. Assume = L 2 (R. 3 ,€ 
Let be a > and A > 0. TTie function & A belongs to L 1 

Proof. As before, we consider as a function of |fc| to simplify the notation. 
We shall prove an estimate of the form 



\B\(r)\ < - 



5^ 



< 



+ r' e 
C A log(2 + r) 

I 1 r 2+2c 



< 



I _)_ r l+2e ' 



< 



1 1 r 2+2e 



(95) 



(96) 



where e > is given by (I13| . The result will follow using a formula similar 
to (l92|) . We first notice that v < w, so that 7?(w + w) — r;(u> — u) > as ( is 
nonincreasing. Hence, by (|13|) . 



1 + f > 1 



rj(w) 



> c 



Inserting in we obtain 

\Bi(r)\<C f - 
Jo t 1 



r 2 1 

T + *2 



+ 1 



t 2 « 



il-2e(l^ + l) J 



< 



c 



1 + r 



2c ' 



(97) 



(98) 



For the three first derivatives of B^, we invoke the following 



Lemma 19. We have for any p — 1,2,3 



1 



d" 



drP \1 + *(r,t,u) 



< 



C 



(1 + *(r, i, u)) (1 + (1 - u)pvp) 



(99) 



Assuming Lemma IT9"! holds, we can write 



(B{)'(r) = -(27T)" 



1 - 1^ 



(l + r 2 t 2 /4) 2 i l + *(r,*,ti) 



fit 



hence 

|(si)'(r)|<c 



o t(l + r 2 t 2 /A) 
1 ri 1+2e di 



/" 2 d { 1 

X (1_u) ^(t+¥m7 

1 e^dt f 1 (l-u)du 



du. 



< 



(r^ + l) 
C 



2+e 



c 



l+e 



1 + (1 - u)r 



2 _|_ r l+2e ■ 

The proof of (j9"6")l is similar. Therefore, we omit it. Instead we turn to the 
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Proof of Lemma \19[ We have 
d ( 1 



B?*(r,t,«) 



dr \ 1 + *(r, t, it) 7 (1 + *(r, i, u)) 2 : 
so that we have to prove that 

C 



1 + *(r,t,u) 



< 



1 + (1 - u)r* 



Since \I/(r, i, u) = ^(il( w + v) + r)(w — v)) + w'u(rj(w + v) — rj(w — v)), we obtain 

\S)(r t lA = 

dr 



i 

\t (r, f, u) = — { (w'+v')r]' (w+v) + (w f — v')rf (w—v)}+w' / u{'i](w+v)—ri(w—v)} 



+ w'u{(w' + v')r)'(w + v) - O' - v')rii'{w - v)} . (100) 

Next we remark that v' = u < 1, w" = 0, to' < 1/2, w" < (7(1 + r)" 1 and 
to'" < C(l + r)~ 2 . Using 1 + W(r, t, u) > 1 + ?7(w + i;) + 77(1/; — u), and Assumption 
(fBf , we obtain an estimate of the form 



< C 



1 



1 



1 



1 + r 1 + 1 10 + u| 1 + \w - v\J ' 



1 + *(r,t, u) 
Eventually, we use 

r t , 

I to + t>| = u> + v > w > — , and \w — v\ = w — v > — (1 — u). 

The proof for the other derivatives is similar. 



We end this section with 

Proposition 20. Assume $j = L 2 (R 3 ,C 4 ) and that ( ^ satisfies ([T2"])-(fT 

We ftaue, as A — > 00 , 



□ 
□ 



b a(°) = ^ log A + 0(1). 



2 

3^ 

Proof. Taking z = \/l — i 2 , we obtain 



z 2 -z 4 /3 



Wo (i-- 2 )(i + c(x^j)) 



dz. 



As C > 0, we have 



1 f E(A) 



z 2 - z 4 /3 



(1-z 2 ) 1 + C 



A a (l-z a ) 



dz < - 

7T 



1 /"■but Z 2 - Z 4 /3 



1-Z 2 



dz = B° A {0). 



To get a lower bound, we use that £ is smooth, and write that ((x) < cx for 
any < x < 1 and some c > 0. We obtain 



1 r e(a) 



z 4 /3 



Wo (i-- 2 )( 4 + c(a^))) 



dz > 



1 /" E(A) 



z 4 /3 



dz 



> 



TT^l - C/A 2 7o 



Ax/l-c/A 2 

B(A) z 2 _ z 4/ 3 log(l + c) 

^-dz = B° (A) - 6V _ + o(l), 



1-z 2 



3tt 
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as ((x) < cx for any < x < 1 and some c > 0. Finally, by (fT3|) , 



(z 2 - z 4 /3)dz 



< CE(AY 



dz 



(l_ a a)i- e /a 



0(1). 



^(i-* 2 )(i + c(a^)) 

This yields the result. □ 
A.3 Study of b A for ((t) = t 

We finally turn to the special cut-off (,(t) = t which was used in the study of 
the ionization in Theorem [U Formula (|94|) yields in this case 



B{(k) = it- 1 



dt 



t{l + \k\H*/A)J 



1-u 2 



l+ir( 4 



|k|" 



3|fe| 2 



(101) 



Notice that B A is nonnegative, when A > l/y/2, and 



4(0) = tt 1 



1 -u 2 



-C?7i = 7T 



1 ^2 



z 2 (l -z 2 /3)rfz 

A 2 -l ~ 
A 2 Z 



A 2 (l-z 2 ) 

A 2 (3A(2A 2 - 3)arctanh (^^) + (8 - 5A 2 )VA 2 - l) 
~ 9%/A 2 - 1(A 4 - 2A 2 + 1) 

= A log A - — + A log 2 + 0(A~ 2 log A). 
dir \)ir Sir 

Hence 5^(0) = B A (0) + 0(A~ 2 log A). Moreover, it can be seen that B A (0) < 
2/(3n) log A when A > 4. The main result of this section is 

Proposition 21. Assume ft = i 2 (R 3 ,C 4 ) and ((t) = t 2 . Let be a > and 
A > 1 such that aB A (0) < 1. The function b A satisfies 

aB{(0) 



< 



LHm l-aB A {0) 
Proof. It follows from (TTUTj) that 

nl 8A 2 dt f VT ^ J (1 - u 2 )du 1 



(102) 



B{(k) 



TTt 3 



l + 3u 2 Ml (i) 2 + |fc| 2 ^ 2 (t,u) 2 + \k\ 2 

where m(t) = 2/t, and fi 2 (t,u) = 2A(1 - 1/A 2 + l/t 2 ) 1 / 2 (l + 3m 2 )- 1 / 2 . The 
Fourier inverse of (a* 2 + |fc| 2 ) -1 is the Yukawa potential e'^ /(4tt\x\) > 0. 
Therefore, the Fourier inverse /a = J-^ 1 {B A ) is nonnegative, so that 



f A (x)dx= / \f A (x)\dx = (2Tr)-V 2 B{(0). 



(103) 



In particular, the operator T : g s L 1 (M. 3 ) i— * T 1 (aB A g) is bounded by \\T\\ < 
<xB A (0). Hence, 1+T is invertible when aB A (0) < 1, and 



ki + TT 1 !! < 



l-afii(0)' 



Proposition |2"T1 follows using that i> A = T(l + T) _1 . 



□ 
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